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We derive the formalism for the leading order corrections to the adiabatic approximation to the scattering 
of composite projectiles. Assuming a two-body projectile of core plus loosely-bound valence particle and a 
model (the core recoil model) in which the interaction of the valence particle and the target can be neglected, 
we derive the non-adiabatic correction terms both exactly, using a partial wave analysis, and using the eikonal 
approximation. Along with the expected energy dependence of the corrections, there is also a strong dependence 
on the valence-to-core mass ratio and on the strength of the imaginary potential for the core-target interaction, 
which relates to absorption of the core in its scattering by the target. The strength and diffuseness of the core- 
target potential also determine the size of the corrections. The first order non-adiabatic corrections were found 
to be smaller than qualitative estimates would expect. The large absorption associated with the core-target 
interaction in such halo nuclei as 11 Be kills off most of the non-adiabatic corrections. We give an improved 
estimate for the range of validity of the adiabatic approximation when the valence-target interaction is neglected, 
which includes the effect of core absorption. Some consideration was given to the validity of the eikonal 
approximation in our calculations. 
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I. INTRODUCTION 

Halo nuclei are often described by few-body models which 
recognise selected internal degrees of freedom of the nucleus. 
Few-body reaction models assume the projectile nucleus con- 
sists of core and valence clusters, which interact with the tar- 
get while also interacting with each other. Using this few- 
body approach, excitations of the projectile to both bound and 
continuum states may be included within the reaction model; 
this breakup is significant for weakly bound systems, such as 
halo nuclei. Further approximations are usually necessary to 
handle the breakup continuum. 

One such approach which includes strong couplings to the 
continuum is the adiabatic approximation [|IJ], in which the 
breakup continuum is collapsed onto a single channel which 
is degenerate with the ground state of the projectile. This is 
often referred to in the Coulomb excitation literature as the 
sudden approximation. Another approach to solving the few- 
body problem, without requiring the use of the adiabatic ap- 
proximation, is the Coupled Discretised Continuum Channels 
(CDCC) method. Here, convergence of the observables may 
be obtained in most cases without the need for further approx- 
imations. CDCC methods, in which the breakup continuum is 
truncated and discretised, have been effectively applied to re- 
actions involving two-body projectiles, such as deuterons and 
single neutron halo nuclei. However, the adiabatic approach 
provides a simpler calculation scheme and as yet, is the only 
method for dealing with multi-channel effects, including the 
continuum, for three-body projectiles. Three-body models are 
frequently used for two-neutron halos such as 6 He, 11 Li, and 
14 Be. 
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A special case of the adiabatic model is when the interac- 
tion between the core and the target dominates and the in- 
teraction between the valence particle and the target can be 
neglected [^]; we call this the core recoil model, since the 
only way the projectile can be excited is through recoil of the 
core in its scattering by the target. This has been applied to 
the elastic scattering of halo nuclei, such as n Be, where the 
core-to-valence mass ratio is large [||]. The core recoil model 
is also particularly applicable to Coulomb dominated process 
when the valence particle is neutral, and has been applied to 
the Coulomb breakup of high energy deuterons [j^, | ] and for 
one and two neutron halos [pi ||| . 

With the current interest in radioactive nuclear beams op- 
erating at energies in the region of tens of MeVs, the valid- 
ity of high energy models, such as the adiabatic model, at 
lower energies is of importance. Comparisons of CDCC and 
adiabatic methods for deuteron induced reactions have been 
extensively studied P]. The adiabatic model was shown to 
be accurate at much lower energies than one would expect 
from qualitative estimates based on the basic assumption of 
the adiabatic approximation — that the breakup energies of the 
projectile excited during the reaction are small relative to the 
centre of mass energy of the projectile. A simple estimate of 
the accuracy of the adiabatic approximation for elastic scat- 
tering and elastic breakup is given in Ref. [^J within the core 
recoil model. However, this estimate did not take into con- 
sideration processes such as core absorption, which has been 
suggested as an important factor in improving the accuracy of 
the adiabatic approximation [Q]. 

The adiabatic approximation underlies many of the micro- 
scopic theories for nuclear reaction, such as Glauber theory 
[B]. Glauber models have been extensively applied to reac- 
tions involving deuterons and light halo nuclei [|l0| [H]] and 
have been useful for extracting information on halo sizes due 
to their diffuse nature [12 |T^]. Much work has gone into 
correcting the eikonal assumptions made in Glauber theory 
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Jl4| , O]. Corrected calculations compare well to few -body 
calculations which make only the adiabatic approximation, 
based on the method of Refs. ^ ^ . 

The formulae for the leading order non-adiabatic correc- 
tions in the core recoil model were given first in Ref. fll^]. 
This paper derived the leading order corrections using the 
eikonal approximation, and discussed the relevance of the 
terms involved. Here, we calculate these corrections for 
n Be and 6 He elastic scattering from a 12 C target at 10 
MeV/nucleon. The non-adiabatic corrections were found to 
be identically zero for a pure point Coulomb potential JI^]. 
The cross terms for the case of nuclear+Coulomb interactions 
are also small for light ion reactions [p0|]. The accuracy of the 
adiabatic approximation for small angle Coulomb scattering 
was discussed in Ref. [||. In this paper, all Coulomb effects 
are ignored. It is felt that this provides an adequate framework 
for discussing the validity of the adiabatic approximation, al- 
though no comparison with experimental data is possible. A 
full comparison of the adiabatic approximation with experi- 
mental data, for the elastic scattering of a single neutron halo, 
including Coulomb effects, has been made in Ref. [ph. 

The format of the paper is as follows. In Section II we 
discuss the adiabatic approximation in the special case of the 
core recoil model, where the valence-target interaction is ne- 
glected. The validity of the adiabatic approximation in the 
core recoil model is examined using qualitative arguments, 
and comparisons with CDCC calculations are made. In Sec- 
tion |l| we rederive the formulae in Ref. using an alter- 
native method which expresses the internal Hamiltonian op- 
erator for the projectile as the product of two operators. We 
provide an eikonal derivation as in Ref. [g|] and then go on to 
derive the corrections exactly using a partial wave sum. The 
eikonal approximation provides useful insights into the nature 
of the non-adiabatic corrections while the exact calculation is 
used to examine the range of validity of the eikonal approxi- 
mation to the corrections. This provides some justification for 
the use of the eikonal approximation in more complete cal- 
culations which include the valence-target interaction. This 
will be reported elsewhere [59, ^jh. In Section |vj we present 
calculations of the non-adiabatic correction in the core recoil 
model for n Be and 6 He. Section gives an improved esti- 
mate for the range of validity for the adiabatic approximation. 



A summary and concluding remarks are given in Section VI 



II. ADIABATIC APPROXIMATION 

Here we describe the projectile using a two-body model. 
The position of the valence particle (v) relative to the core (c) 
is described by the vector r. The position vector of the centre 
of mass of the whole projectile from the target is R. The few- 
body Hamiltonian is 



H = T R + V cT (R - or) + V vT (R ~ M + H v 



where 



T R = h 2 K 2 R /2n \K R = 
U 2 K 2 r /2ii vc + V vc {r) 



(1) 



(2) 
(3) 



Here, H vc is the internal Hamiltonian of the projectile which 
has a ground state wavefunction, <fio, satisfying 



(H vc + e )</»o = 0, 



(4) 



where £o is the binding energy of the projectile. The core- 
target and valence-target potentials are V c t and V v t respec- 
tively, while V vc is the core-valence potential that binds the 
projectile. The reduced masses of the projectile-target and 
core-valence systems are /i and [i vc respectively, while the 
mass ratios a = m v /mp and (3 — —m c /mp = a — 1 relate 
to the distance of the core and valence clusters from the center 
of mass of the projectile. 

The exact wavefunction for the few-body scattering system 
is [63 



I* 



(+)\ 

K I 

E+ 



E+ - H 
E + ie 



(5) 
(6) 



where E is the total energy of the few-body scattering system 
and K is the centre of mass momentum. The adiabatic ap- 
proximation ^ assumes that the internal motion of the pro- 
jectile is frozen, and thus the internal Hamiltonian of the pro- 
jectile can be replaced by a constant. This is chosen to be the 
ground state energy of the projectile (— £o)> because then the 
incident part of the three-body wavefunction has unit ampli- 
tude. The adiabatic wavefunction is therefore 



I* 



ad(H 
K 



Eft -T R - V cT - V vT 



\4>o,K), (7) 



where Eq = E + eq is the centre of mass energy of the pro- 
jectile. 



A. Core recoil model 

A special case of the adiabatic model is when the scatter- 
ing is dominated by the core-target interaction. This is the 
core recoil model in which the valence-target interaction is 
neglected. This provides a considerable simplification of the 
adiabatic Hamiltonian, as the only place that the vector r now 
appears is in the core-target potential. This dependence can 
be transformed away using the translation operator 



Ur(x) 



_ -ix-K R 



(8) 



The wavefunction in the core recoil model is 



^>(R,r) = <t> (r)e >ar -"x ( K J (R-<xr), (9) 



ar-K(+) , 



where Xk '' is me two-body distorted wave for a particle of 
mass /i in the core-target potential, 

(R\ X £ ) ) = {R\ieG ( £\K), (10) 

and is the adiabatic Green's operator with the valence- 
target potential switched off, and contains the core-target po- 
tential in the R co-ordinate only: 



1 



£ + -T R - V c t(R) 



(11) 
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The core recoil model wavefunction, when evaluated in the 
appropriate few-body T-matrix with the valence-target inter- 
action neglected, allows a factorisation into a T-matrix for a 
point particle of mass n, multiplied by a formfactor, Foo M\: 

T$(K,K') = {K'\V cT \x ( k ] )F^Q) (12) 
F m {aQ) = (0 o |e iQT - Q |0o>. (13) 

The elastic differential cross section can then be obtained from 
that of a point particle multiplied by the formfactor squared: 

The point particle elastic cross section is that obtained for a 
particle of reduced mass // interacting via the core-target po- 
tential, with a centre of mass energy Eq. The formfactor in- 
cludes all effects of excitation and breakup of the halo struc- 
ture. 



TABLE I: Woods-Saxon potential parameters for core-target and 
core-valence interactions. Energies are in MeV and lengths in fm. 
The range parameters are to be multiplied by the appropriate masses. 



c+T 


V 


Rv 


ay 


W 


Rw aw 


10 Be+ 12 C 
q+ 12 C 


123.00 
37.16 


0.750 
1.846 


0.800 
0.452 


65.00 
13.27 


0.780 0.800 
1.846 0.452 


c+v 


V 


Rv 


av 






10 Be+n 
a+2n 
a+2n 


86.42 
172.17 
134.82 


1.000 
0.800 
0.800 


0.530 
0.300 
0.728 


(ground state) 
(d-wave resonance) 



Using values for the range and diffuseness of the core-target 
interaction from Table |[ some typical values for A are 

n Bc + 12 C at 10 MeV/nucleon => A = 0.14, 
6 Hc + 12 C at 10 MeV/nucleon => A = 5. 



B. Validity of the adiabatic approximation 

The range of validity of the adiabatic approximation has of- 
ten been investigated by comparison with CDCC calculations 
which do not make the adiabatic approximation, but instead, 
use a discretised continuum spectrum of breakup energies for 
the projectile [j7|. By examining the relative excitation and 
breakup energies required to gain convergence in the CDCC 
calculations it is found that the adiabatic approximation gives 
a better estimate of the elastic cross section than one might 
expect. 

Estimates of the average excitation and breakup energies in- 
volved in the scattering of the projectile have been formulated 
using the core recoil model [||. The adiabatic approximation 
is expected to be valid when the collision time of the projectile 
and target, t co \\, is short in comparison to the time associated 
with the core- valence internal motion, t[ nt . In Ref. [||], an 
upper limit on the ratio of these times was given as 



£coll 



< A, 



where 



A = 2 



Ro 1 



m c (mp + wit) a 2 K 



(15) 



(16) 



and Rq defines the range of the core-target interaction and a 
the diffuseness. The adiabatic approximation is then expected 
to be valid when 



A«l. 



(17) 



Note that A is strongly dependent on the valence-to-core mass 
ratio, and only weakly dependent on the incident energy of 
the projectile, through the 1/K dependence. Also note that 
the derivation of A in Ref. [^J involves an estimate of the exci- 
tation energy involved, but the projectile binding energy does 
not appear. 



These values set upper limits on the time ratio that must be 
much less than unity. From these values we would expect that 
when the valence-target interaction is neglected, the adiabatic 
approximation would give a reasonable description of the ex- 
act cross section for n Be, but a poor description for 6 He, at 
10 MeV/nucleon. 



C. Adiabatic versus CDCC numerical calculations 

Extensive studies of deuteron elastic scattering have shown 
that there is more to the adiabatic approximation than the basic 
assumption that the excitation energy of the projectile is small 
in comparison to the centre of mass energy of the projectile. 

Fig. [j] compares adiabatic and CDCC calculations for 
n Be and 6 He elastic scattering from a 12 C target at 10 
MeV/nucleon. The calculations neglect the valence-target in- 
teraction and do not include any Coulomb potential. The 
adiabatic cross section with the valence-target interaction ne- 
glected is just the core recoil model cross section obtained 
from Eq. Q). 

The core-target potential parameters are given in Table [j| 
The 10 Be radius parameters are to be multiplied by 10 1 / 3 + 
12 1 / 3 and the a by 12 1 / 3 . 

The ground state wavefunction of 11 Be is assumed to be 
a pure 2si/2 neutron single particle state, with a separation 
energy of 0.503 MeV, calculated in a central Woods-Saxon 
potential (Table The potential depth is adjusted to obtain 
the required binding energy of the 10 Be+n system. Assuming 
a core root mean squared (rms) radius of 2.28 fm, this gener- 
ates the n Be composite nucleus, with rms radius of 2.90 fm, 
in agreement with a recent few-body analysis of halo sizes 

@& 

The CDCC calculations were performed using FRESCO 
fl23|]. The lpi/2 bound state in 11 Be, with an excitation en- 
ergy of 320 keV, was included. The 10 Be+n continuum was 
modelled using 10 bins from 0-20 MeV for each of the s-,p- 
,d-, and /-wave breakup states. 
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FIG. 1 : CDCC (solid) versus adiabatic (dashed) calculations for Be 
and 6 He elastic scattering from a 12 C target at 10 MeV/nucleon. The 
valence-target interaction has been neglected as well as the Coulomb 
interaction. The cross section in the limit of no excitation or breakup 
is represented by the dotted line. The 11 Be cross sections are multi- 
plied by 100. 



A two-body di-neutron model was assumed for the ground 
state wavefunction of 6 He, with a binding energy of 0.975 
MeV. The di-neutron is assumed to be in a 2s single parti- 
cle state. The a+2n potential used to generate the two-body 
ground state wavefunction was assumed to have a Woods- 
Saxon form (Table |). 

The a+2n continuum, due to the breakup of 6 He, was mod- 
elled using 15 bins from 0-30 MeV for each of the s-,p-, and 
/-wave breakup states, using the same potential as that used 
for the ground state wavefunction. The ri-wave breakup states 
included a resonance at 1.8 MeV above the ground state, with 
a bin below the resonance and 15 bins up to 30 MeV above 
the resonance. A potential which reproduced a 1.8 MeV res- 
onance in a+2n with a width of 1 13 keV, was used for these 
ri-wave breakup states (Table ||). 

The cross section in the limit of no excitation or breakup is 
calculated by folding the core-target potential over the ground 
state wavefunction, to form a two-body projectile-target in- 
teraction. This is shown by the dotted line as a reference to 
highlight the importance of excitation and breakup channels 
using the different methods. 

Fig. [j] shows that the effect of the adiabatic assumption is 
negligible for 11 Be a nd small for 6 He, even though the esti- 
mate in Section |IIB| suggests otherwise. To understand why 
this is the case, the leading order corrections to the adiabatic 
approximation are evaluated in the following sections. 



III. FIRST ORDER NON-ADIABATIC CORRECTIONS 



tile, for which (H vc + sq) is non-zero. Therefore, we expand 
the wavefunction of Eq. (|5|) in powers of (H vc + Eq), where to 
first order we have 



I* 



(.+)) 

K I 



I* 



ad(+) x 
K I 



G { +\H vc + e Q )\^>). (18) 



ad(+). 



The first order non-adiabatic correction to the elastic adiabatic 
T-matrix is then 



AT$(K,K') = (*p->\(H vc + e )\*™ 
The adiabatic wavefunction is 



ad(H 



where V'jc i s me distorted wave for the two-body 
projectile-target scattering system in the potential V c t(R — 
ar) + V v t{R — (3r), with a reduced mass \i, and normalised 
so that it has an incident plane wave in R with unit amplitude. 

The operator product of the 3 factors 4>o(r), (H vc + Eq), 
and 4>o (r) which appears in Eq. ( |T9| ) can be expressed as (see 
Appendix Q) 



ad(+) 
K 



(R,r) = Mr)iPK (+ \R,r), 



(19) 



(20) 



4>o(r)(H vc + £ )<M r ) 



(21) 



This result assumes that 0o is an s-state and hence, without 
loss of generality, can be assumed to be real. The result also 
assumes that V vc is local. 

The expression ( pT| ) can be conveniently evaluated by al- 
lowing the two operators to operate in opposite directions, i.e., 
the left del on the bra and the right del on the ket. Thus, the 
matrix element is reduced to 



m=±i,o (22) 



A. Non-adiabatic corrections to the core recoil model 

Using the analytical wavefunction of the core recoil model, 
the matrix element for the first order correction to the adi- 
abatic elastic T-matrix can be simplified significantly. By 
translating the two-body distorted waves, 



,(-> 



UrX { k\R) 



X { J(R-ar) 



K 

x ( K kR) 



u 



R- 



(23) 
(24) 



the adiabatic wavefunctions in the core recoil model can be 
written as 



^ d(+) 



(R,r) 



«■ ' • \k (25) 



iar-(K R -K) (+), 

X ( J(i?)Ve iQr < K n~ K '\ (26) 



Corrections to the adiabatic approximation arise because 
the scattering process mixes in excited states of the projec- 



The derivatives with respect to r are now straight-forward 
since the r dependence has been separated. Therefore, when 
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we take the derivative followed by the inner product of tions. Using the eikonal approximation to the distorted waves, 
Eqs. ( j25| ) and (^), the Kr operators in the exponents can- 
cel, allowing the formfactor of Eq. (|l3]) to be factored out. 
This leaves us with the first order correction to the adiabatic 
elastic T-matrix: 



AT; 



ad 



-X(K,K')F 00 (aQ) 



where 



X(K,K') = {xU\(K R - K') ■ (K R - K)\ x £>). (28) 



(+)\ 



(27) 



X ( k\R) = e iK - R e-^l-°° dziVMb+Zl *\ (29) 
e -i«:'-R e -^/~d4VcT(b'+4*') )( 30) 



and by making the small angle approximation, K ■ K' = 1, 
which is consistent with the eikonal approximation, Eq. ( p7| ) 
can be written as Jl9| ] 



B. Evaluation of the non-adiabatic corrections using the 
eikonal approximation 

In Ref. [p9[], Eq. ( p7| ) was evaluated using the eikonal ap- 
proximation. This leads to simple formulae for the correc- 



AT a d = ( X { K ) |AV cT | x W)Foo(«Q), (3D 



where 



J 



d Zl V b V cT (b + Zl K) 
I 



dz 2 V b V cT (b + z 2 K) , 



(32) 



and 7 is the mass ratio 



rriT 



(33) 

m c [mp + niT) 

In evaluating the matrix element in Eq. (f3l|), the eikonal 
distorted waves combine to form the eikonal phase factor, 



XcT(b) = - — 



dzV cT (b + zK), 



(34) 



in which the z integration has been performed. Note that this 
is possible since, for elastic scattering, K' = K. The z inte- 
gration in Eq. (|T]) therefore only involves AV c t, and we can 
define the quantity 



K f°° 
X(b) = -— / dzAV C T- 

^0 J-oo 



(35) 



For a central potential, the integral over all z of AV c t can 
be simplified to a single integral over all z of the potential 
squared, 



dzV? T { \/b 2 + z 2 
We can now write Eq. (^TJ) as a scattering amplitude, 
A/, 



/eik 
00 



Kl bdbJ (Qb)S(b) x (b)F 00 (aQ). 
o 



(36) 



(37) 



The eikonal scattering amplitude in the core recoil model, 
including first order non-adiabatic corrections, is then 

POO 

=iK bdb J (Qb) [1 - S(b)(l + ix)} F 0Q (aQ), 
Jo 

(38) 



where ffg< = / ^ k + A/ e ^. Here, 



S(b) = exp 



-— / dzv cT {^T7 2 ) 

^0 J-oo 



(39) 



is the eikonal S'-matrix for 11 Be as a point particle moving in 
the core-target potential. 



We see that the non-adiabatic correction term, x of Eq. (|3f 
is a dimensionless quantity, whose size relative to unity deter- 
mines the magnitude of the corrections. 

The structure of the expression (^6|) for the correction factor 
X can be traced to the physical origin of corrections to the adi- 
abatic approximation. In the few-body model used here these 
corrections arise from excitation of the projectile through tidal 
forces generated by the interaction of the core with the target. 
These forces arise because of the displacement of the core 
from the centre-of-mass of the projectile and it is therefore 
natural that the mass ratio m v /m c and derivatives of the core- 
target potential should be crucial. There is a quadratic depen- 
dence on V c t because a two-step process must be involved 
if the projectile is to end up in the ground state (elastic scat- 
tering). The extra 1/Eq factor in x, over and above that ex- 
pected from the energy dependence of the eikonal phase (see 
Eqs. ( |34| ) and ([35])), must be related to the ex pecte d depen- 
dence on the collision time discussed in Section II B . 



We note that x is multiplied by the eikonal S'-matrix which 
will restrict the impact parameters contributing to the overall 
corrections. 
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C. Exact evaluation of non-adiabatic corrections 



with asymptotic form (ignoring the Coulomb interaction) 



While the eikonal approximation provides useful insights 
into the nature of the non-adiabatic corrections, the matrix el- 
ement in Eq. ( |28| ) can be evaluated exactly using a partial wave 
expansion. The quantity X(-K", K'), defined in Eq. (^8|), was 
shown in Ref. Jl^ ] to have the alternative form, 



X(K,K') 



-(Xk'\ [Kr, Vct] GlV 2 [K R , V cT ) \ X { k ] ) 



We can write Xk as a P ar tial wave sum, 

,(+) 



(40) 



R\X { K ] ) =4*^i t Y; m (K)Y em (R) X f\R), (41) 



I'm 



where Xz is me solution of the radial Schrodinger equation 

(E o -T e -V cT )xi +) (R)=0, (42) 
and T( is radial kinetic energy operator, 



T, 



ti z 1 d 2 



h 2 £{£+!) 



2(i RdR 2 2fi R 2 

,(+) 



(43) 



ki 1 



(-) 



(44) 



where Si is the partial wave S'-matrix, and are radial 

Hankel functions, i.e., irregular solutions of the equation 



1(1+1) 2(i 



Vi = 0. 



(45) 



dR 2 R 2 

With a similar expression for x^> and using a partial wave 
sum for ' , the matrix element of Eq. ( |40| ) cai 

X(K,K') = 47r^P £ (cos0)X £ , 
where we have defined 



(46) 



i? 2 di? X 



(+) 

if,«-l 



(47) 



The non-adiabatic correction to the elastic scattering ampli- 
tude in the core recoil model is then 



A/oo = -7^oo (qQ) J2 Pt(«*9)Xt 



(48) 



(50) 



and where can take two values, which leads to the two terms 
in Eq. @: 

t' = £±l. (51) 

We have defined Sw as the coefficient of the Hankel func- 
tion in the inhomogeneous solution. It can be shown that it 
is equivalent to the difference between the S'-matrix for the 
homogeneous solutions for I and I', 



Sw — ±{Si — Sf) 
-(+) 



£' = e±i. 



(52) 



The oscillatory nature of X^A, (R), for large R, means that 
the integrals in Eq. (^) have to be dealt with carefully. The 
asymptotic form for the solutions to the differential Eqs. ( |42| ) 
and (49) is reached when R is outside the range of the poten- 



tial. The radial Hankel functions, H^, have the asymptotic 
form 



H, 



(±) 



R^oa 



The partial distorted wave, x\ , has the asymptotic form (ig- when 
noring the Coulomb interaction) 



KR->t(i + \). 



(53) 



(54) 



The solutions of the differential equations were computed out 
to a radius R = R , which is outside the range of the potential. 
Then, the integrals from i?o to Re, where Re is chosen so that 
the asymptotic condition (|54j) is met, are performed using the 
explicit form of the Hankel function [E4J], 



H, 



Z^2>s\{l-a)\ 



(KRY 



(55) 



The integral of the asymptotic form of the Hankel function 
from Re — > oo can be performed analytically by making the 
substitution R = iy + Re, and integrating over the complex 
plane, 



/ &Re 2{( - KR -^^ = i / dye- 2K y+ 2i 

JR e , JO 



KR.i-it'ir 



1 2iKR lt -ii'w 
2K 



(56) 



It it well known that a partial wave sum can be written ex- 
actly as an integral over impact parameters |^5j|. By making 
the semi-classical correspondence, £ = bK, and assuming the 
scattering angle is small, the sum over partial waves can be 
written as [ j2o] l 

y* oo 

VP £ (cos6») -> K I dbJ (Qb). (57) 
i Jo 



In Eq. (pj[), Xkm' is me solution of the inhomogeneous 



equation 



{Eo-T if -V^r)X^,(R) 



<49» 



The eikonal equivalent of Xt (Eq. 
substituting this relation into Eq. 
Eq. (pi) to obtain 



)) can then be found by 
and comparing it with 



jXt = bS(b)x(b), 



(58) 
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where 7 is the mass ratio of Eq. (33), which also appears on 
the R.H.S. within the correction term, \. It is then under- 
stood that it is the magnitude of either side of Eq. ( p8| ) which 
determines the overall size of the non-adiabatic corrections: 
it contains the overlap of the correction term with the elastic 
scattering S-matrix, which determines which impact parame- 
ters contribute to the cross section. 



IV. NUMERICAL RESULTS 

A. Application to 11 Be+ 12 C elastic scattering 

Here we examine the non-adiabatic corrections in the core 
recoil model for n Be scattering from a 12 C target. n Be is a 
good example of a single neutron halo nucleus to which the 
core recoil model can be applied. The interaction between 
the 10 Be core and the target dominates the elastic scattering 
and it is thus a reasonable first approximation to neglect the 
neutron-target interaction. The core recoil model was applied 
to the elastic scattering of 11 Be+ 12 C at 49.3 MeV/nucleon in 
Ref. [0]. At this energy we expect the adiabatic approxima- 
tion to be good, and we therefore calculate the non-adiabatic 
corrections at 10 MeV/nucleon, where we have calculated the 
qualitative estimates for the validity of the approximation ear- 
lier. 



T he p otentials and wavefunctions were discussed in Sec- 



tion 



1 II C . The potential used for the 10 Be core in Ref. [|zp was 
obtained from elastic scattering data at 59.4 MeV/nucleon. 
Due to the absence of data at lower energies, the same po- 
tential will be used at 10 MeV/nucleon, but as the correc- 
tions are dependent on the potential geometry, ideally the 
potential should be fixed by elastic scattering of 10 Be at 10 
MeV/nucleon. 

In Fig. U (bottom), the correction term, x, is plotted against 
impact parameter for 11 Be+ 12 C at 10 MeV/nucleon. It has a 
maximum value of approximately 0.6, which is a significant 
correction as it is added to unity, but it is then multiplied by the 
eikonal S'-matrix, which is plotted on the same impact param- 
eter scale in Fig. || (top). We see that the non-adiabatic correc- 
tions are largest for small impact parameters, but the S'-matrix 
is zero in this region, due to a large imaginary component in 
the core-target potential. This kills off most of the correc- 
tion term, with only large impact parameters, corresponding 
to grazing collisions, contributing to the cross section. The 
large impact parameters correspond to forward scattering an- 
gles where there is little momentum transferred to the projec- 
tile during the scattering, and therefore only small corrections 
to the adiabatic approximation. This can be seen more clearly 
in Fig. ||, where the magnitude of the S-matrix and correction 
term, x, is plotted along with the overlap of the two functions 
(multiplied by 10). It is the maximum size of this overlap, 
|Sx|max, in comparison to unity, which determines the over- 
all size of the non-adiabatic corrections. We see that for this 
system, the maximum value is 0.015, so the non-adiabatic cor- 
rections are small. 
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FIG. 2: Eikonal S'-matrix and non-adiabatic corrections for 
11 Be+ 12 C in the core recoil model. In the top figure, the solid line 
represents the two-body S-matrix for 1 Be as a point particle inter- 
acting with the 12 C target via the core-target interaction. The dotted 
line includes the first order non-adiabatic corrections through multi- 
plying the S-matrix by the factor (1 + ix), while the correction term, 
X, is plotted in the bottom figure on the same impact parameter scale. 



B. Accuracy of eikonal calculations 

The eikonal approximation provides us with an understand- 
ing to the nature of the non-adiabatic corrections, but the use 
of the eikonal approximation must be validated, as in the en- 
ergy region we are considering, we would expect significant 
non-eikonal corrections. 

The non-a diabatic corrections were formulated exactly in 
Section III C ; by comparison of the exact non-adiabatic cor- 
rections with those calculated in the eikonal approximation, 
the validity of the eikonal approximation for these calcula- 
tions can be assessed. 

Eq. ( |58| ) shows that jXt can be compared to the overlap of 
the correction term and the S-matrix in the eikonal approxi- 
mation, as shown in Fig. ||, multiplied by b. 

Each side of Eq. ( jBlf ) is plotted in Fig. [|, with the exact cal- 
culation (L.H.S. of Eq. (58)) represented by the circles, and 

T8t)) represented by the 



the eikonal calculation (R.H.S. of Eq. 
line. The exact calculation is plotted for each partial wave and 
scaled to match the corresponding impact parameter for the 
eikonal calculations. This comparison shows that the overlap 
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FIG. 3: Overlap of S'-matrix and non-adiabatic corrections for 
11 Be+ 12 C. The dashed line is the magnitude of the correction term \ 
plotted against impact parameter. The dot-dashed line is the magni- 
tude of the .S'-matrix. The solid line is the magnitude of the product 
of the S'-matrix and the correction term x multiplied by 10. 
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FIG. 4: Exact and eikonal calculations of non-adiabatic correc- 
tions for 11 Be+ 12 C at 10 MeV/nucleon. The solid line represents 
the eikonal approximation to the non-adiabatic corrections plotted 
against impact parameter along the bottom axis. It is the overlap of 
the correction term and the S'-matrix weighted by the corresponding 
impact parameter. The circles represent the exact non-adiabatic cor- 
rections for each partial wave, which are labelled along the top axis 
and scaled to match the corresponding impact parameter via the re- 
lation I — bK, where K = 3.97 fm _1 . The exact correction term, 
Xi, is scaled by the mass ratio 7. 



FIG. 5: Angular distribution of the elastic differential cross section 
for n Be+ 12 C in the core recoil model with no Coulomb interac- 
tion. The upper figure compares the cross sections in the core recoil 
model, using exact (solid) and eikonal (dashed) calculations of the 
two-body cross section. The lower figure shows the non-adiabatic 
correction as a fraction to the cross section, with the lines corre- 
sponding to exact (solid) and eikonal (dashed) calculations. 



of the eikonal S'-matrix and the correction term, x> gives a rea- 
sonable representation of the non-adiabatic corrections. The 
corrections are slightly over-estimated by the eikonal approx- 
imation, especially for the smaller impact parameters, but the 
larger impact parameters are well repreduced. This is where 
we expect the eikonal approximation to do better as this cor- 
responds to smaller scattering angles. 

In Fig. U (top), the differential cross section in the core re- 
coil model is plotted. The two curves represent the different 
methods of calculating the two-body cross section in the core 
recoil model. The solid line uses an exact partial wave analy- 
sis while the dashed line uses the eikonal approximation; the 
formfactor is the same in both cases. In the bottom of Fig. || 
the first order non-adiabatic correction is plotted as a fraction 
to the core recoil cross section. We see that the fractional cor- 
rection to the adiabatic approximation is well reproduced by 
the eikonal approximation, even though the eikonal cross sec- 
tion differs significantly from the exact cross section in the 
core recoil model. The eikonal approximation over-estimates 
the cross section at this energy, but we saw from Fig. [| that 
it also over-estimates the magnitude of the corrections. As 
a fraction to the cross section, however, the non-adiabatic 
corrections are well reproduced. The large core absorption, 



9 



which kills of most of the corrections, means that only large 
impact parameters contribute to the cross section, and there- 
fore, the eikonal approximation gives a reasonable description 
of the non-adiabatic corrections. 

We note that the fractional corrections to the adiabatic ap- 
proximation do not depend on the formfactor, and therefore 
do not depend on the internal structure of the projectile. This 
is a result of our assumption of an s-wave projectile which re- 
sults in the factorisation shown in Eq. (^7). For a non-s-wave 
projectile, the additional terms in Eq. (Al) contribute. 

an estimate on the accuracy of the adia- 



1.2 



In Section II B 



batic approximation was given as an upper limit on a time 
ratio, which had to be much less than one. This upper limit 
was calculated to be 0.14 for n Be+ 12 C at 10 MeV/nucleon: 
we then would expect the adiabatic approximation to be good 
at this energy from this estimate. From the calculations of 
the first order corrections, the adiabatic approximation is ex- 
tremely accurate for thi s sys tem, going beyond the range that 
the estimate of Section II B, due to the key role of the strong 
absorption associated with the scattering at small impact pa- 
rameters. The maximum overlap function had a value of 0.015 
at the peak, which suggests that the adiabatic approximation 
is approximately 10 times better that the qualitative estimate 
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FIG. 6: Overlap of .S'-matrix and non-adiabatic corrections for 
6 He+ 12 C. The dashed line is the magnitude of the correction term 
X plotted against impact parameter. The dot-dashed line is the mag- 
nitude of the S'-matrix. The solid line is the magnitude of the product 
of the S-matrix and the correction term \- 



C. Application to He+ C elastic scattering 

In the previous section, the 11 Be+ 12 C system was studied 
because it was a reasonable approximation to use the core 
recoil model, as the ratio of the valence-to-core masses was 
1/10. This small ratio also meant that the corrections in this 
model were small. The large absorption in the 10 Be+ 12 C po- 
tential also played an important role in accuracy of the adia- 
batic approximation for elastic scattering. To examine the role 
of the mass ratio and core absorption, the core recoil model 
was applied to 6 He+ 12 C elastic scattering. 

The 6 He nucleus has a two-neutron halo with an a core, 
so the ratio of valence-to-core masses is 1/2. The a core is 
light and will appear slightly transparent to the 12 C target, so 
corrections at small impact parameters will contribute. Elas- 
tic scattering cross sections are available over a wide range 
of energies for a+ 12 C, so a more realistic potential can be 
used. The potentials and wavefunctions were discussed in 
Section II C . 

The magnitude of the non-adiabatic correction term, x, is 
plotted in Fig. ^ (dashed line). We see that the 6 = value is 
approximately the same as for 11 Be at the same energy. De- 
spite the valence-to-core mass ratio being five times larger for 
6 He, the smaller a+ 12 C potential cancels this out, leaving a 
correction term similar to that for 11 Be. The magnitude of 
the peak at the potential surface for x is larger than its b = 
value; for n Be it was smaller. This is because the a+ 12 C po- 
tential has a sharper potential surface, increasing the deriva- 
tive of the potential at the surface, thus increasing the correc- 
tion term (see Eq. (|36|)). By comparison with the S-matrix 
(dot-dashed line), we see that the peak in the correction term 
now appears in a region of the S-matrix which is non-zero; 



therefore, the overlap with the S'-matrix is more significant. It 
has a maximum overlap of around 0.3, compared to 0.015 in 
the 11 Be case. This factor of 20 difference between the two 
cases arises even though the magnitude of x at b = is of the 
same order in both cases. 

The overlap between the correction term, x, and the S'- 
matrix is increased in the 6 He case because the core-target 
potential has a weaker absorption associated with it. This in- 
creases the magnitude of the S-matrix in the region of max- 
imum corrections, thus increasing the overlap. Even so, the 
S-matrix is still relative small at the peak of the correction 
term, x, and so the maximum overlap is still much smaller 
than the maximum size of the correction term. The maximum 
overlap of 0.3 is small in comparison to unity, thus corrections 
to the adiabatic approximation are still small, even though the 
qualitative estimate in Section II B suggested that the approx- 
imation would be poor at this energy. 

The accuracy of eikonal calculations for the 6 He case is 
shown in Fig. [7| We see, as in the 11 Be case, that the non- 
adiabatic corrections are over-estimated by the eikonal ap- 
proximation, with the worst agreement for smaller impact pa- 
rameters, but we see in Fig. |] that this difference is more ev- 
ident in the fractional correction to the cross section. The 
S-matrix is not as strongly absorbing as for n Be, so these 
small impact parameters are contributing to the cross section 
at large scattering angles. For the largest scattering angles 
(above 25°), the eikonal approximation to the first order non- 
adiabatic corrections reduce the cross section, while an exact 
evaluation shows an increase in the cross section. The eikonal 
approximation was good for 11 Be because the large core ab- 
sorption meant that only large impact parameters contributed. 



10 



10 



He+ C E h =60MeV 



15 



20 



25 

yX, (exact) 
bSx (eikonal) 




FIG. 7: Exact and eikonal calculations of non-adiabatic correc- 
tions for 11 Be+ 12 C at 10 MeV/nucleon. The solid line represents 
the eikonal approximation to the non-adiabatic corrections plotted 
against impact parameter along the bottom axis. It is the overlap of 
the correction term and the S'-matrix weighted by the corresponding 
impact parameter. The circles represent the exact non-adiabatic cor- 
rections for each partial wave, which are labelled along the top axis 
and scaled to match the corresponding impact parameter via the re- 
lation i — bK, where K = 2.77 fm _1 . The exact correction term, 
Xi, is scaled by the mass ratio 7. 



This is not the case for 6 He scattering and so the eikonal ap- 
proximation is not as accurate. 
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FIG. 8: Angular distribution of the elastic differential cross section 
for 6 He+ 12 C in the core recoil model with no Coulomb interac- 
tion. The upper figure compares the cross sections in the core recoil 
model, using exact (solid) and eikonal (dashed) calculations of the 
two-body cross section. The lower figure shows the non-adiabatic 
correction as a fraction to the cross section, with the lines corre- 
sponding to exact (solid) and eikonal (dashed) calculations. 



D. Dependence on core absorption 

In the previous section, we see that the non-adiabatic cor- 
rections 6 He+ 12 C at 10 MeV/nucleon are larger than for 
11 Be+ 12 C at the same energy. The reasons for this are two 
fold: firstly, the core-to-valence mass ratio is much larger for 
the former; secondly, the core-target potential has a smaller 
imaginary component for the former, due to less absorption of 
the core. The second effect is examined in more detail here. 

The potential we have used for the 10 Be+ 12 C interac- 
tion was obtained at 59.4 MeV/nucleon. As the non- 
adiabatic corrections have been calculated at the energy of 10 
MeV/nucleon, we would expect that the imaginary potential 
strength to be reduced, but without any experimental data at 
this energy, its precise value cannot be fixed. 

The dependence of the corrections on the imaginary po- 
tential strength, Wo, is shown in Fig. ^[ The non-adiabatic 
corrections are plotted for a range of imaginary potential 
strengths for the 11 Be+ 12 C reaction at 10 MeV/nucleon. 
Fig. |9| contains three graphs: the top graph plots the modulus 
of the correction term, x; the middle figure is the modulus of 
the S'-matrix; the bottom figure is the overlap of the two. The 



dependence of the corrections on the imaginary potential is 
shown for various imaginary potential strengths: the solid line 
represents the 65 MeV imaginary potential that was obtained 
from the elastic scattering of 10 Be+ 12 C at 59.4 MeV/nucleon, 
the dashed line corresponds to a 40 MeV imaginary potential, 
the dotted line is for Wn=20 MeV, and the dot-dashed line 
represents Wn=10 MeV 

We see that different imaginary potential strengths do not 
effect the correction term \ significantly, but the effect on the 
S'-matrix is large. As the imaginary potential strength is re- 
duced, the correction term, x, is reduced slightly; but, the tar- 
get appears more transparent and so the S-matrix is increased 
significantly for the small impact parameters. The overlap of 
the correction term with the S-matrix is then significantly in- 
creased, as shown in the bottom figure, and it is this overlap 
which determines the overall size of the corrections. 

The accuracy of the eikonal approximation for evaluating 
the non-adiabatic corrections was shown to be suspect for 
6 He, which was said to be due to the weak imaginary poten- 
tial. To see this effect for the 11 Be case, the 20 MeV imag- 
inary potential was used to compare the exact and eikonal 
calculations of the non-adiabatic corrections. The fractional 
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FIG. 9: Non-adiabatic correction term for 11 Be+ 12 C at 10 
MeV/nucleon for varying imaginary potential depths. The top figure 
plots the modulus of the correction term, x, versus impact parameter 
while the middle figure shows the modulus of the S'-matrix on the 
same impact parameter scale. The bottom figure is the overlap of the 
S'-matrix and \- The imaginary potential depths for the core-target 
potential are given in the legend. 



non-adiabatic corrections are shown in Fig. |l0| We see here 
that the eikonal approximation reproduces the exact correc- 
tions poorly for large scattering angles (above 20°), as in the 
6 He case. The eikonal approximation fails to predict the in- 
crease in the cross section for large scattering angles that was 
seen when the exact non-adiabatic corrections were included. 
The non-adiabatic corrections for n Be are still small com- 
pared to the 6 He calculations, even when the imaginary poten- 
tial depth is reduced, due to the smaller valence-to-core mass 
ratio. 



V. IMPROVED ESTIMATE FOR VALIDITY OF 
ADIABATIC APPROXIMATION 

We have shown in our numerical evaluation of the non- 
adiabatic corrections that strong core absorption improves the 
accuracy of the adiabatic approximation, going beyond that 



expected from the estimates of Section II B 



The key point is that when strong absorption is present, it 
is the value of the correction term, \ (Eq. d3~6[)), in the re- 
gion near the strong absorption radius which is of importance. 
The size of the non-adiabatic corrections is determined by the 
maximum overlap of \ ar, d the S'-matrix. If you model % by 
an exponential and the S'-matrix by a Woods-Saxon, then a 
good estimate for the size of the non-adiabatic corrections is 
X, evaluated at the strong absorption radius, multiplied by one 
half (the value of the S-matrix at the strong absorption radius). 

We can estimate the value of \ (Eq. ( 36 )) at the strong ab- 
sorption radius as the potential in this region has the form of 
a simple exponential, so that 



dzV 2 T {Vb 2 + z 2 ) = / 

-OO J — ( 



dzV^e 



2 -2Vb 2 +z' 2 /a 



(59) 



where a is the diffuseness of the core-target interaction. Since 
the most important values in the z integral are those around 
z = 0, the square root can be expanded in powers of z/b: 
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FIG. 10: Angular distribution of the elastic differential cross section 
for n Be+ 12 C in the core recoil model with the imaginary potential 
depth reduced to 20 MeV. The non-adiabatic corrections are plotted 
as a fraction to the cross section, with the lines corresponding to exact 
(solid) and eikonal (dashed) calculations. 



dz y2 e -2Vb^/a 



&zVo 2 e -2i,(1+ * /2b )/o - 



(60) 

The integral of the square of the potential can therefore be 
written 



dzV? T ( \Jb 2 + z 2 ) « Vb^V 2 T (b), (61) 
and an estimate of x, near the strong absorption radius, is 

X(b) « 2^^±b" 2 V 2 T {b). (62) 



db 



The derivative of the potential dominates in -^b 3 / 2 V^ T {b), 
and as the potential only depends on b, the overlap of the cor- 
rection term and the S-matrix at the strong absorption radius 
can be written, 



\Sx\i 



WE 2 



(63) 



R=R B 
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where R s is the strong absorption radius, and 7 is the mass 
ratio in Eq. (33). The adiabatic approximation, for scattering 
systems with strong core absorption, is then valid when 



\Sx\ t 



< 1. 



(64) 



For strong absorbing systems, this criterion replaces condi- 
tion (|F|h 

Eq. (|53[) gives approximate values for |Sx| max of 0.01 and 
0.2 for n Be and 6 He respectively. This compares well with 
calculated values of 0.015 and 0.3 from Figs. [}] and |[ T hese 
values greatly improve on the estimates given in Section |1I B| 
and give the maximum overlap to the correct order of magni- 
tude for scattering systems with core absorption. 



VI. SUMMARY 



to be of great importance in the contribution from these differ- 
ent processes. The corrections for 11 Be, although being very 
small partly due to the small valence-to-core mass ratio, could 
be much larger when valence particle recoil is included. This 
will be dealt with elsewhere pflKlJ], 
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We have calculated first order non-adiabatic corrections 
for the first time using the core recoil model, in which the 
valence-target interaction is neglected. Two reactions were 
studied: the elastic scattering of n Be and 6 He from a 12 C tar- 
get at 10 MeV/nucleon. The non-adiabatic corrections were 
compared to previous qualitative estimates of the validity of 
the adiabatic approximation. 

The eikonal approximation was used to gain insights into 
the nature of the non-adiabatic corrections. They were found 
to be dependent on the overlap of a correction term (Eq. (36)) 
and the S-matrix. The correction term was found to be 
strongly dependent on the ratio of the valence mass to that 
of the core, as with the qualitative estimates. Along with the 
expected energy dependence, there was also a dependence on 
the strength and diffuseness of the core-target interaction. The 
overlap of the correction term with the S'-matrix was strongly 
dependent on the strength of the imaginary potential for the 
core-target interaction. Strong core absorption kills of most of 
the non-adiabatic corrections as the maximum of the correc- 
tion term lies in a region where the S'-matrix is zero, whilst 
smaller imaginary potentials increase the overlap of the of 
the S-matrix with the correction term producing larger non- 
adiabatic corrections. The corrections calculated are much 
smaller than what is expected from qualitative estimates due 
to the key role that core absorption plays. 

An improved estimate for the validity of the adiabatic ap- 
proximation, when the valence-target interaction is neglected, 
is given in Section |v| which includes the role of core absorp- 
tion. This new estimate recognises that when the core-target 
absorption is strong, it is the size of the correction term in the 
region near the strong absorption radius which is of impor- 
tance in determining the size of the non-adiabatic corrections. 
Eq. ( |63| ) gives a value which is of the correct order of mag- 
nitude for the size of the non-adiabatic corrections when core 
absorption is present. 

In the core recoil model we have used here, non-adiabatic 
corrections can only arise through projectile excitations occur- 
ring from recoil of the core in its scattering by the target. If 
the valence-target interaction is to included, corrections could 
also arise through recoil of the valence particle. The strong 
dependence on the valence-to-core mass ratio is then expected 



APPENDIX A: OPERATOR IDENTITY 

We are interested in the operator product of the three factors 
ijj*(r), (H — e), and <f>{r). Our result is that this can be re- 
expressed as 



r(H ~ e)4>{r) = - %-V r tl>*</> ■ V r 



(Al) 



+ 2~ [(VW*)4>-V>* (V r </>)]• V r . 

The result assumes that <j> is an eigenstate of H = — j- V 2 + V 
with eigenvalue e, but ip can be arbitrary. V must be a local 
operator. <fi and ip do not have to be bound states. 

For an s-wave state <j> in a real potential, tj> can be assumed 
to be real and the last term in Eq. ( Al ) vanishes if%j) = <p: 



(r)(H - e)0(r) 



2ft 



(A2) 



In proving Eq. ( \1 ) we use round brackets to indicate when 
the V operator acts only locally on the functions inside the 
brackets. We have 

V r VW r = (V r yj*<jy) ■ V r + (^*0)V 2 

+ V*V, 2 - ^*(V 2 0) - 2^*(V r 0) • v r 
- [(V r ^*)0-^*(V r 0)] • v r 

+ ^*V 2 0-^((^ -£)</>), (A3) 



where we have put double brackets around the last factor to 
emphasise that V acts on <fi only. We can drop these brackets 
if V is a local operator. 
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The last line in Eq. ( 43 ) can now be re-expressed in terms 
of H and we obtain 



where 



V^-Vr = [(V P ^>-V*(V P 0)]-V P 



V-e] 



ft 



This is equivalent to the identity given in Eq. ( Al ). 



(A4) 



The identity (A.2) is frequently seen in a form which is 
equivalent to it when <p(r) is nodeless: 



AW = 



-1 



(A6) 
(A7) 



H-e = — 4>- l v r 4? ■ v r 4>~ 1 

2/x 



2/2 



(A5) 



and p = — iftV,. is the momentum operator. The factorised 
form ( A5) is used in 1 -dimension in the formulation of super- 
symmetric quantum mechanics 
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